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Abstract 

The Los Alamos version of the Quark-Gluon String Model implemented in the code LAQGSM03 was extended 
to describe photonuclear reactions at energies up to tens of GeV. We have incorporated 56 channels to consider 
jp elementary interactions during the cascade stage of reactions and 56 channels for 771 interactions, in addition 
to absorption of photons on two-nucleon pairs, into LAQGSM03 and have tested the extended code referred to as 
LAQGSM03.01 on a number of measured high-energy photonuclear reactions. Here, we present the photonuclear 
reaction model of LAQGSM03.01 and show examples of its results compared with available experimental data. 
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1. Introduction 



Following an increased interest in nuclear data for such projects as the Accelerator Transmutation of 
nuclear Wastes (ATW), Accelerator Production of Tritium (APT), Spallation Neutron Source (SNS), Rare 
Isotope Accelerator (RIA), Proton Radiography (PRAD) as a radiographic probe for the Advanced Hydro-test 
Facility and others, for several years the US Department of Energy has supported our work on the development 
■ of improved versions of the Cascade- Exciton Model (CEM) of nuclear reactions QlEl and the Los Alamos version 
of the Quark-Gluon String Model (LAQGSM) [3J 0] to describe reactions induced by particles and nuclei at 
energies up to about 1 TeV/nucleon jEj-[H]. To describe fission and production of light fragments heavier 
than 4 He, we have merged both our codes with several evaporation/fission/fragmentation models, including the 
Generalized Evaporation/fission Model code GEM2 by Furihata ^Ojj the fission-like binary-decay code GEMINI 
by Charity et al. and the Statistical Multifragmentation Model (SMM) by Botvina et al. Our codes 

perform as well as and often better than other current models in describing a large variety of spallation, fission, 
and fragmentation reactions, therefore they are used as event-generators in several transport codes. 

However, the initial version of the Los Alamos Quark-Gluon String Model (LAQGSM) just like its 
precursor, the Quark-Gluon String Model (QGSM) did not consider photonuclear reactions, while its 2003 
version, LAQGSM03 pp, describes such reactions only for energies up to about 1.5 GeV. This is not convenient 
when using our codes to solve problems for PRAD, NASA, and other high-energy applications or to analyze 
future high-energy measurements at the Thomas Jefferson National Accelerator Facility (CEBAF), where pho- 
tons with much higher energy are created and need to be simulated with our even-generators in transport 
codes like MCNP6 [H], MARS HS|, MCNPX [THj, and others. To address this problem, we have extended 
[S] LAQGSM03.01 to describe photonuclear reactions up to tens of GeV. The next section presents a short 
description of the high-energy photonuclear reaction model of LAQGSM03.01, followed by two sections with 
examples of calculated particle spectra, yields, and recoil properties of products from high-energy photonuclear 
reactions compared with available data. 

2. High-Energy Photonuclear Model of LAQGSM03.01 

Initially , an improved version of the Dubna InraNuclear cascade photonuclear reaction model developed 
originally 37 years ago by one of us in collaboration with Iljinov and Toneev |17j to describe photonuclear 
reactions at energies above the Giant Dipole Resonance (GDR) region was incorporated into LAQGSM. [At 
photon energies T 7 = 10-40 MeV, the de Broglie wavelength X is of the order of 20-5 fm, greater than the average 
inter-nuclconic distance in the nucleus; the photons interact with the nuclear dipole resonance as a whole, thus 
the INC is not applicable.] Below the pion production threshold, the Dubna INC considers absorption of photons 
on "quasi-deuteron" pairs according to the Levinger model |18|: 

T Z{A-Z) 

0-jA = L — CTyd , (dj 
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where A and Z are the mass and charge numbers of the nucleus, L ss 10, and cr 7 d is the total photoabsorption 
cross section on deuterons as defined from experimental data. 

At photon energies above the pion-production threshold, the Dubna INC considers production of only one 
or two pions (channels #1-3 and 15-17 from Table 1); the concrete mode of the reaction is chosen by the Monte 
Carlo method according to the partial cross sections, defined from available experimental data. This limits the 
use of such an approach to photon energies of only up to about 1.5 GeV, as contribution from multiple pion 
production become predominant at higher energies. 

To address this problem, we have extended LAQGSM03 pQ to describe photonuclear reactions at energies 
up to tens of GeV. For this, we took advantage of the high-energy event generators for ^yp and 777 elementary 
interactions from the Moscow high-energy photonuclear reaction model 19 kindly sent us by one of its coauthors, 
Dr. Igor Pshenichnov. Let us note that the 7p and 771 event generators from the Moscow INC ^5] use 
an improvement and extension to higher energies of the Genova jp and 771 event generators developed by 
P. Corvisiero et al. to describe photon- nucleon interactions up to 4 GeV [201 ■ We have incorporated into 
LAQGSM03.01 56 channels to consider jp elementary interactions during the cascade stage of reactions, and 56 
channels for 771 interactions. These reaction channels new to LAQGSM03.01 are listed in Table 1 (its precursor, 
LAQGSM03 pQ, considers only channels #1-3 and 15-17, and only up to about 1.5 GeV). 

To describe in LAQGSM03.01 the two-body channels #1-14, we use part of a file containing smooth ap- 
proximations through presently available experimental data sent us by Dr. Pshenichnov but have developed 
our own algorithms and written our own routines to simulate unambiguously da/dQ and to choose the corre- 
sponding value of O for any E 11 using a single random number £ uniformly distributed in the interval [0,1], as 
described in pQ . Fig. 1 shows examples of angular distributions of 7r° from ^p — > 7r°p interactions as functions 
of 0^ s at eight different photon energies as simulated by our routines compared with available experimental 
data. Examples of similar distributions for tt + from jp — > ir + n interactions may be found in pp. 

To describe the channels #15-21 with two and three pions in the final state, we use the jp and 777 event 
generators send us by Dr. Pshenichnov, but use our own interpolation for integral cross sections. We do not 
show here examples of cross sections for these channels, as we reproduce with LAQGSM03.01 all the results by 
Pshenichnov et al. shown in Figs. 6 and 8 of Ref. |19j . 

Finally, to describe in LAQGSM03.01 the multi-pion channels #22-56, we use the isospin statistical model 
as realized in the 7p and 777 event generators send us by Dr. Pshenichnov and described in details in |19| . 
without any changes. For channels #22-56, we reproduce exactly in LAQGSM03.01 the results by Pshenichnov 
et al. as published in Ref. |19|. therefore we do not show here examples of such results. 

After the bombarding photon is absorbed by two nucleons or interacts inelastically with a nucleon according 
the channels #1-56, we get inside the nucleus several "secondary" cascade nucleons, pions, or other mesons and 
resonances listed in Table 1, depending on which channel is simulated from the corresponding cross sections at 
the given photon energy to actually occur. These "secondary" cascade particles interact further with inranuclear 
nucleons or leave the nucleus, depending on their coordinates and momenta. The further behavior of the reaction 
starting from this stage, after the photon had "disappeared", is described by LAQGSM03.01 exactly the same 
way as for any other types of reactions, induced, e.g., by nucleons or heavy ions. 

LAQGSM03.01 [H] is the latest modification of the Los Alamos version of the Quark-Gluon String Model 
p5], which in its turn is an improvement of the Quark-Gluon String Model |13| . It describes reactions induced 
by both particles and nuclei as a three-stage process: IntraNuclear Cascade (INC), followed by preequilibrium 
emission of particles during the equilibration of the excited residual nuclei formed after the INC, followed by 
evaporation of particles from or fission of the compound nuclei. The INC stage of reactions is described with 
a recently improved version of the time-depending intranuclear cascade model developed initially at Dubna, 
often referred in the literature simply as the Dubna intranuclear Cascade Model, DCM (see j2J an d references 
therein). The preequilibrium part of reactions is described with the last version of the Modified Exciton Model 
(MEM) from the improved Cascade-Exciton Model (CEM) code CEM03.01 0. The evaporation and fission 
stages of reactions are calculated with an updated and improved version of the Generalized Evaporation Model 
code GEM2 by Furihata ^Ojj which considers evaporation of up to 66 types of different particles and light 
fragments (up to 28 Mg). If the excited residual nucleus produced after the INC has a mass number A < 13, 
LAQGSM03.01 uses a recently updated and improved version of the Fermi Break-up model (in comparison with 
the initial version developed in the group of Prof. V. S. Barashcnkov at JINR, Dubna and used in QGSM J3| 
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lels of elementary 7iV interactions taken into account in LAQGSM03.01 
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and described in 22 ) to calculate its decay instead of considering a preequilibrium stage followed by evaporation 
from compound nuclei, as described above. LAQGSM03.01 considers also coalescence of complex particles 
up to 4 He from energetic nucleons emitted during the INC, using an updated coalescence model in comparison 
with the version described in [21]. 

We have developed also "S" and "G" versions of LAQGSM03.01, namely: LAQGSM03.S1 and LAQGSM03.G1 
0. LAQGSM03.S1 |Hj is exactly the same as LAQGSM03.01, but considers also multifragmentation of excited 
nuclei produced after the preequilibrium stage of reactions, when their excitation energy is above 2A MeV, 
using the Statistical Multifragmentation Model (SMM) by Botvina el al. ^2] (the 'S" in the extension of 
LAQGSM03.S1 stands for SMM). LAQGSM03.G1 |U is exactly the same as LAQGSM03.01, but uses the 
fission-like binary-decay model GEMINI of Charity et al. ^3 > which considers evaporation of all possible frag- 
ments, instead of using the GEM2 model [TUj (the "G" stands for GEMINI). 
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Figure 1. Example of eight angular distributions of ir° from 7p — > ir°p as functions of m s at photon energies 
from 260 MeV to 1.4 GeV. The dashed lines show the old approximations of the Dubna INC 17 while the 
solid lines are our new approximations incorporated into LAQGSM03.01 (and into CEM03 and LAQGSM03 
PQ). References to experimental data shown by different symbols may be found in our recent work pQ. 
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3. Particle Spectra 



We have tested LAQGSM03.01 and its "S" and "G" versions extended to describe high-energy photonuclear 
reactions against practically all experimental data above 1 GeV we were able to find in the literature, but we 
limit ourselves here to presenting only several exemplary results. 

Let us start with comparing results by LAQGSM03.01 with experimental data at relatively low energies, of 
about 1 GeV, where other photonuclear models also do work, then move gradually to reactions at higher energies. 
Figure 2 presents examples of proton spectra from bremsstrahlung interaction with carbon at E max = 1050 
MeV. One can see that LAQGSM03.01 describes as well as CEM03.01 does the measured proton spectra and 
agrees with the data better than the direct knockout model and the quasi-deuteron calculations [35] do. 

Figure 2. Comparison of measured |23| differential cross 
section for proton photoproduction on carbon at 43°, 90°, 
and 154° by bremsstrahlung photons with E max = 1.05 
GeV (symbols) with CEM03.01 (solid histograms) and 
LAQGSM03.01 (dashed histograms), and predictions by 
the direct knockout model |24| (dotted lines) and a quasi- 
deuteron calculation |25| (dot-dashed lines), respectively. 
The experimental data and results by the direct knockout 
and quasi-deuteron models are taken from Fig. 5 of Ref. |24| . 

Figure 3 presents examples of proton and pion spectra 
from photonuclear reactions at higher energies, namely 
proton spectra at 60, 90, and 150 degrees from interaction 
of bremsstrahlung 7 quanta of maximum energy E max = 
2.0, 3.0, and 4.5 GeV with 12 C, 27 Al, 63 Cu, and 208 Pb and 
pion spectra from interaction of the same bremsstrahlung 
7 quanta with 12 C, 63 Cu, and 208 Pb. These reactions 
are of a higher interest to us than the one presented in 
Fig. 2 for the following reason: The proton and pion 
spectra presented in Fig. 3 were measured more than 25 
years ago at the Yerevan Physics Institute [23-123 with 
a major hope of revealing some "exotic" or unknown 

mechanisms of nuclear reactions leading to the production of so called "cumulative" (i.e., kinematically forbidden 
for quasi- free intranuclear projectile- nucleon collisions) particles. At the time the measurements [2I]-!2I1 were 
done, there were not high-energy photonuclear reaction models available in the literature, therefore these data 
were not analyzed so far by any models. We do not know of any publication or oral presentation where these 
measurements were reproduced by a theoretical model, event generator, or transport code. It is noteworthy that 
LAQGSM03.01 describes quite well all the spectra measured both in the cumulative and non-cumulative regions 
(as it does, e.g., for particles emitted from interaction of 400 GeV protons with nuclei; see for details) in a 
single approach, without any fitting or free parameters, and without involving any "exotic" reaction mechanisms. 

These results do not imply, of course, that the 7- or proton-nucleus interaction physics is completely 
described by the reaction mechanisms considered by LAQGSM03.01. Our present results do not exclude 
some contribution to the production of these particles from other reaction mechanisms not considered by 
LAQGSM03.01. But the contribution from "exotic" mechanisms to cumulative particle production from these 
high-energy reactions seems to be small; inclusive particle spectra are not sensitive enough for an unambiguous 
determination of the mechanisms of particle production, just as observed heretofore at intermediate and low 
energies |3l7) . 

4. Product Yields and Recoil Properties 

This section presents several examples of cross sections (yields) of the nuclides produced in high- 
energy photonuclear reactions and of their recoil properties. Fig. 4 shows a comparison of the measured 
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Figure 3. Proton spectra at 60, 90, and 150 degrees from interaction of bremsstrahhmg 7 quanta of maximum 
energy E max = 2.0, 3.0, and 4.5 GeV with 12 C, 27 Al, 63 Cu, and 208 Pb (left and right top panels and left bottom 
panel of four plots). Spectra of 7r + and ir~ produced by E max = 4.5 GeV bremsstrahlung on 12 C and spectra 
of charged pions (both 7r + and ir~) from interaction of the same bremsstrahlung 7 quanta with 12 C, 63 Cu, 
and 208 Pb (right bottom panel of four plots). Experimental values shown by symbols are from j2E]-|2nj while 
histograms show results by LAQGSM03.01. To the best of our knowledge, we are able to describe these data 
with LAQGSM03.01 for the first time (see text). 
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Figure 4. Detailed comparison be- 
tween experimental yields |31| and those 
calculated by LAQGSM03.01 of radioac- 
tive products from the interaction of 
bremsstrahlung 7 quanta of maximum en- 
ergy 4.5 GeV with 93 Nb. The cumulative 
yields are labeled as "c" while the inde- 
pendent cross sections, as "i" . To the best 
of our knowledge, using LAQGSM03.01 
we are able to describe these data for the 
first time. 



production cross sections from in- 
teraction of bremsstrahlung 7 quanta 
of maximum energy 4.5 GeV with 93 Nb 
EH with results by LAQGS03.01. We 
see that LAQGSM03.01 describes most 
of the measured yields quite well, with 
only several cases when the calculations 
differ from the data within a factor of 
two to three, that is also not too bad, 
considering that all calculations are in 
the absolute value, without any fitting 
of parameters or normalization. 
Fig. 5 shows a similar reaction, induced by E max — 4.5 GeV bremsstrahlung on 65 Cu, but measured 

regarding the recoil properties of the produced nuclides |32] . No product yields were measured from this 
reaction, so the top-left plot in Fig. 5 shows only predictions by LAQGSM03.01 and its "S" and "G" versions 
for this quantity. The left plot in the middle row of Fig. 5 shows predictions by LAQGSM03.01 and its 
"S" and "G" versions for the mean laboratory angle of products. We see quite a big difference for this 
characteristics between predictions by versions of the model with different treatment of the evaporation stage 
of reaction (compare the "standard", LAQGSM03.01, results that uses GEM2 with results by its "G" version 
that uses GEMINI) and between results calculated without (LAQGSM0.01) and with taking into account 
multifragmentation of highly-excited compound nuclei (LAQGSM03.S1). Unfortunately, this quantity was not 
measured and we can not uncover the "real" reaction mechanisms based on these results until experimental 
data are available. 

The other three plots in Fig. 5 show results by LAQGSM03.01 and its "G" and "S" versions for the Z- 
averaged A-dependence of the mean laboratory velocity v z of products, the R=F/B ratio of the forward product 
cross sections to the backward ones, and their mean laboratory kinetic energy compared with available data 
|32| . We see a very good agreement between calculations and the data for the mean kinetic energy, a not 
so good agreement for the mean parallel velocity v z of products, and a poor agreement for R=F/B. Let us 
note that in the experimental work |32|. these characteristics were measured only for several final isotopes, 
while results of calculations shown in Fig. 5 are for the Z-averaged A-dependence of all possible products. 
If to compare the data with calculations for only the measured products, isotope-by-isotope, the agreement 
between measurements and calculations becomes better. However, as can be seen from Fig. 6, there is still some 
disagreement between experimental data and theoretical results by all three versions of LAQGSM considered 
here for the values of R=F/B, similar to results obtained for reactions induced by protons and deuterons |33| . 
In making such comparisons, we first recognize that the experiment and the calculations differ in that: 1) the 
experimental data were extracted assuming the "two-step vector model" (see references and details in |32II33| L 
while the LAQGSM calculations were done without the assumptions of this model; 2) the measurements were 
performed on foils (thick targets), while the calculations were done for interactions of photons (protons and 
deuterons, in the case of Ref. [221) with nuclei (thin targets). These differences must be considered before 
assessing possible deficiencies in the models. 
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Figure 5. Results by 
LAQGSM03.01 and its "S" 
and "G" versions for the 
Z-averaged A-dependence 
of the mass yield of all 
products, their mean par- 
allel laboratory velocity v z 
in the beam direction, their 
mean laboratory angle O, 
the R=F/B ratio of the 
forward product cross sec- 
tions to the backward ones, 
and their mean laboratory 
kinetic energy for the re- 
action of E max = 4.5 GeV 
bremsstrahlung photons 
on 65 Cu compared with 
available experimental data 
by Arakelyan el al. |32) . 

Finally, Fig. 7 
shows a comparison 
of the calculated by 
LAQGSM03.01 mass 
distribution of products 
with experimental data 
for the highest energy we 
were able to find data in 
the literature, namely for 
the reaction E max = 6 
GeV bremsstrahlung 
gammas on 238 U [521155] . 
We see a good agreement 
between the calculations 
and these high-energy 
photonuclear reaction 
data. 



To conclude, let us note that to the best of our knowledge, we were able to describe with LAQGSM03.01 all 
the reactions discussed in this section for the first time: We do not know of any publication or oral presentation 
where these experimental data were reproduced by a theoretical model, event generator, or transport code. 



5. Summary 

The latest modification of the Los Alamos Quark-Gluon String Model released in the code LAQGSM03.01 
and its "S" and "G" versions were extended to describe photonuclear reactions at energies up to tens of GeV. 
We have tested our high-energy photonuclear models against practically all experimental data above 1 GeV 
we were able to find in the literature and can conclude that they describe the measured data reasonably well, 
without re-fitting any parameters. Our models are not yet ready to describe properly reactions induced by 
gammas of tens of MeV, in the GDR region (though they provide not completely wrong results even at such 
low energies; see, e.g., pQ). We are developing now approximations for the GDR region, to extend the use of 
our codes for low-energy photonuclear reactions as well. LAQGSM03.01 was incorporated recently into the 
transport code MARS ^ and is currently being incorporated into MCNP6 [TJ and MCNPX [THj. This would 
allow others to use LAQGSM03.01 as an event-generator in these transport codes to simulate high-energy 
photonuclear reactions with targets of practically arbitrary geometry and nuclide composition. 
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Figure 6. The measured data |32| for the mean laboratory 
kinetic energy and the R=F/B ratio of the forward product 
cross sections to the backward ones for the same reaction as 
shown in Fig. 5, but compared with results by LAQGSM03.01 
and its "S" and "G" versions for only the measured products 
(listed on the top of both plots) , isotope-by-isotope. 
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Figure 7. Comparison of measured (filled circles) 
|34l I35| mass distribution of the nuclides produced by 
E max = 6 GeV bremsstrahlung photons on 238 U with 
results by LAQGSM03.01 (thick solid line). Open circles 
show some additional experimental data measured at 
other energies and compiled from the literature in Tab. 
5 of Ref. EU; results by LAQGSM03.01 at E max = 1.0 
to 7.4 GeV are shown for comparison by thin and thick 
dashed lines, respectively. 
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